Surface Tension Prediction for Liquid Mixtures

Recently Escobedo and Mansoori (AIChE J. 42(5): 1425, 1996) proposed a new expression, originally derived from statistical mechanics, for the surface tension prediction of pure fluids. In this report, the same theory is extended to the case of mixtures of organic liquids with the following expression, σm = [(1-Tr,m) Tr,m exp (0.30066/Tr,m + 0.86442 Tr,m) (P 0,m ρL,m P 0,m ρv,m)] 4 P0,m = { ΣiΣj [xixj (Pcij / Tcij) 7/3 P 0,ij]} 1/4 [ΣiΣj (xixj Tcij / Pcij)] 7/12 P0,ij=(1-mij)(P0,iP0,j) where σm is the surface tension of the mixture. Tr,m=T/Tc,m and Tc,m is the pseudo-critical temperature of the liquid mixture. ρ L,m and ρ v,m are the equilibrium densities of liquid and vapor, respectively. P0,m and P V 0,m are the temperature-independent compound-dependent parameters for the liquid and vapor, respectively. Pc,i , Tc,i and xi are the critical pressure, critical temperature and the mole fraction for component i, respectively. P0,i is the temperatureindependent parameter for component i and mij is the unlike interaction parameter. Using the proposed equation surface tensions of 55 binary mixtures are predicted within an overall 0.50 AAD% which is better than the other available prediction methods. When the P0,i’s are made compound-insensitive using a corresponding states expression the surface tension of all the 55 binary mixtures are predicted within an overall 2.10 AAD%. The proposed model is also applicable to multicomponent liquid mixtures. ______________________________________________________________________________ (*) To whom correspondence concerning this article should be addressed. AIChE Journal, Vol. 44, No. 10, pp. 2324-2332, 1998 2 Introduction The surface tension of a liquid mixture is not a simple function of the surface tensions of the pure liquids. Also, the composition of the bulk phase and the composition at the vapor-liquid interface are not always the same. At the interface, there is migration of the species having the lowest surface tension, or free energy per unit area, at the temperature of the system. This migration at the interface results in a liquid-phase rich in the component with the highest surface tension and a vapor phase rich in the component with the lowest surface tension (Teixeira et al., 1992, Croxton, 1986, Carey et al, 1980;Winterfield et al., 1978; Massoudi and King, 1975, Wertheim,1976). According to Chapela et al., (1977) the more volatile component is adsorbed from the mixture and extends about one atomic diameter into the gas phase beyond the other species. In most instances the bulk composition is known. However, the composition at the vaporliquid interface is unknown. Therefore, there is a need for relating the surface tension of mixtures to bulk-phase concentrations and properties (i.e. densities). In general, there have been several approaches for estimating the surface tension of mixtures (i) based on empirical or semi-empirical thermodynamic relations suggested first for pure fluids and (ii)based on statistical mechanical grounds. In a recent publication Escobedo and Mansoori (1996) presented a set of new expressions for predicting the surface tension of pure fluids. These new expressions were developed from the statistical mechanical basis of the Macleod equation. It was shown that their new expressions represent rather accurately the surface tension of a variety of pure fluids over a wide range of temperatures. In this paper, the expressions proposed originally by Escobedo and Mansoori (1996) are extended to the case of mixtures. Theory For the surface tension of pure organic fluids, σ, Escobedo and Mansoori (1996), proposed, σ = [P (ρl ρv)] 4 (1) P = P0 (1-Tr) Tr Exp(0.30066/Tr + 0.86442 Tr) (2) where P is known as parachor, ρl is the liquid density, ρv is the vapor density, Tr is the reduced temperature and P0 is a constant which is a function of the characteristics of the molecule under consideration. This new expression, derived from statistical mechanics, represents the experimental surface tension of 94 different organic compounds within 1.05 AAD%. In addition, P0 = 39.6431 [0.22217 2.91042× 10-3 (R/Tbr)] Tc/Pc (3) AIChE Journal, Vol. 44, No. 10, pp. 2324-2332, 1998 3 was proposed as a corresponding-states expression for the temperature-independent, compounddependent parameter P0. In this equation, R * = Rm/Rm,ref; Rm is the molar refraction given by Rm = (4πNA/3)[α +μ 2f(T)] = (1/ρ)(n2-1)/(n2+2); α and μ are the polarizability and the dipole moment of the molecule; NA is the Avogadro’s number, n is the sodium-D light refractive index of a liquid at 20¡ C; Rm,ref=6.987 is the molar refraction of the reference fluid (methane). When this generalized expression is used, surface tensions for all the 94 pure compounds can be predicted within an overall 2.57 AAD% at all temperatures investigated. Exact statistical mechanical extension of equations 1-3 to mixtures results in quite complex expressions. However, using the principles of the conformal solution theory (Massih and Mansoori, 1983) the following generalization of Eq. 2 for the case of mixtures will be concluded, σm = [(1-Trm) Tr,m exp (0.30066/Tr,m + 0.86442 Tr,m) (P 0,m ρL,m P 0,m ρv,m)] 4 (4) where σm is the surface tension of the mixture; P L 0,m and P V 0,m are the temperature-independent conformal parameters for the liquid and vapor phase, respectively; Tr,m is the pseudo-reduced temperature for the mixture (i.e. Tr,m = T / Tc,m); Tc,m is the pseudo-critical temperature of the liquid mixture); and ρL,m and ρv,m are the equilibrium densities of the liquid and vapor phase, respectively. The aim of this report is the application of Eq. 4 for the prediction of the surface tension of mixtures, over the whole concentration range from the temperature-independent parameter P0 for the pure components. In order to apply Eq. 4 to mixtures we need the following, 1. Accurate prediction of phase equilibrium properties (i.e. composition of the liquid and vapor phases, vapor and liquid densities, and saturation pressures). 2. Mixing rules for the temperature-independent parameter P0 in order to calculate the temperature-independent parameters P0,m and P V 0,m for the liquid-and vapor-phase mixtures. Prediction of Phase Equilibrium Properties In this work, the equilibrium properties for a multicomponent mixture were calculated in two steps. First, the Peng-Robinson (PR) equation of state was used to predict the vapor pressure and liquidand vapor-phase compositions. The PR equation of state was selected for this purpose because its parameters are known for many pure compounds and because its interaction parameters are available for numerous binary mixtures. Since the constants of the PR equation are fitted to vapor pressure data it is reliable for vapor-liquid-equilibrium (VLE) calculations. Second, the predicted compositions and vapor pressure by the PR equation were in turn used AIChE Journal, Vol. 44, No. 10, pp. 2324-2332, 1998 4 with the Riazi-Mansoori (RM) equation of state (described below) to predict the densities of the liquid and vapor phases. The RM equation of state yields more accurate density predictions than the PR equation. The calculations described above have been incorporated into a comprehensive computer package for vapor-liquid-equilibrium (VLE) calculations. This computer package has been used to perform the routine bubble-point and flash calculations. The PR equation of state is given by, P = RT / (V b) aα /(V2 + 2bV b2) (5) where P, T, and V=1/ρ are pressure, temperature, and molar volume, respectively; R is the universal gas constant; a = 0.45724 RTc / Pc; b = 0.07780 RTc / Pc; Tc is the critical temperature; Pc is the critical pressure; α = [1 + (0.37464+ 1.5422 ω 0.26992 ω2)(1 Tr )] ; ω is the acentric factor; Tr=T/Tc is the reduced temperature. The mixing rules for the temperature-dependent term (a α) and constant b are given by: (a α)m = ΣΣ [xixj(a α)ij] and bm = Σ (xi bi), respectively; (a α)m is the temperature-dependent parameter for the mixture; bm is the volume-related constant for the mixture; and xi is the mole fraction of species i in the mixture. The cross-parameters are calculated using the following expression, (a α)ij = (1-kij)[(aα)i(aα)j] ; where kij is a binary interaction parameter. This interaction parameter is calculated using a correlation recently proposed by Gao, 1992: (1 kij) = [2 (Tci Tcj) / (Tci +Tcj)] Zcij ; Zcij =[(Zci + Zcj)/2]; where Zci is the critical compressibility factor for pure component i. Bubble-point and flash calculations are performed using the above equation of state following the standard algorithms reported in textbooks (Walas, 1985). The predicted saturation pressure and vaporand liquid-phase compositions are then used along with the Riazi-Mansoori (1993) equation of state (RM EOS) to calculate the equilibrium densities for the vapor and liquid phases. The RM EOS is given below, P = ρRT/(1-bρ) -aρ2/[T1/2.(1+bρ)] (6) a = 0.42748 RTc/Pc and b = (0.08664 RTc/Pc).δ(R ,Tr) δ-1 = 1+{ 0.02[1 0.92 exp (-1000 Tr 1 )]-0.035(Tr 1)} (R* 1) R* = Rm/Rm,ref; Rm = (4πNA/3)[α + μ 2f(T)] = (1/ρ)(n2-1)/(n2+2); Rm,ref=6.987 where ρ is the molar density; T, P, Tc, Pc, Tr, and R are as previously defined; Rm is the molar refraction; Rm,ref is the molar refraction of the reference compound(methane); α and μ are the polarizability and the dipole moment of the molecule; NA is the Avogadro’s number, n is the sodium-D light refractive index of a liquid at 20¡ C (Riazi and Mansoori, 1993). In order to apply Eq. 6 to the case of mixtures the van der Waals mixing rules and the AIChE Journal, Vol. 44, No. 10, pp. 2324-2332, 1998 5 molecular expressions for the cross-parameters (Riazi and Mansoori, 1993) have been used, Tcm = ΣiΣj [xixj (Tcij 2 / Pcij)] / ΣiΣj [xixj (Tcij / Pcij)] (6a) Tcij = (1 kij RM) (Tci Tcj) 1/2 (6b) Pcm = ΣiΣj [xixj (Tcij 2 / Pcij)] / { ΣiΣj [xixj (Tcij / Pcij)]} 2 (6c) Pcij = 8 Tcij / [(Tcii / Pcii) + (Tcjj / Pcjj)] (6d) Rm = ΣiΣj (xixj Rij) (6e) Rij = [R*ii + R*jj] / 8 (6f) In these equations, Tcm is the pseudo-critical temperature of the mixture; Tci is the critical temperature of pure component i in the mixture; Pcm is the pseudo-critical pressure of the mixture; Pci is the critical pressure of pure component i in the mixture; Pcij is a cross-parameter for pressure; kij is the binary interaction parameter (RM stands for Riazi-Mansoori to differentiate it from the interaction parameter used in the PR equation of state); Rm is the dimensionless molar refraction of the mixture; Ri is the dimensionless molar refraction of pure component i in the mixture; and xi is the mole fraction of component i in the mixture. The binary interaction parameter, kij, is calculated using the following expression, (1 kij) = (1 kij) { 8 [(Tcii / Pcii)(Tcjj / Pcjj)] / [(Tcii / Pcii) + (Tcjj / Pcjj)]} (6g) Mixing Rules for the Temperature-independent parameter Po This section is devoted to the development of the mixing rules for the temperatureindependent parameter Po. These mixing rules may be derived considering the proportionality of Po to the intermolecular parameters(or critical properties) according to Kwak and Mansoori (1986). Following the principle of corresponding states a reduced surface tension may be defined as σr = σ / [Tc (kTc)] and expected to be a universal function of Tr = T/Tc (i.e. the reduced temperature). Thus, it can be shown that, Po,ii ∝ Tc 13/12 / Pc (7) It may be assumed that the temperature-independent parameter, Po,m, for the mixture will follow the same proportionality with respect to the critical temperature and pressure of the mixture, Po,m ∝ .Tc,m 13/12 / Pc,m (8) AIChE Journal, Vol. 44, No. 10, pp. 2324-2332, 1998 6 where Tc,m and Pc,m are the pseudo-critical temperature and pressure for the mixture, respectively. Introducing the van der Waals mixing rules for Tc,m and Pc,m, Eqs. 6a and 6c, into Eq. 8 the following expression is obtained, P o,m∝ {ΣiΣj(xixjTcij 2/Pcij)/ΣiΣj(xixjTcij/Pcij)} 13/12/{ΣiΣj(xixjTcij 2/Pcij)/[ΣiΣj(xixjTcij/Pcij)] } (9) this equation may be further reduced to the following expression, Po,m∝ [ΣiΣjxixj(Tcij 2 / Pcij)] [ΣiΣjxixj(Tcij / Pcij)] 7/12 (10) the term (Tcij 2 / Pcij) maybe further simplified to obtain, Tcij 2 / Pcij =(Tcij 11/12 / Pcij ) (Tcij 13/12 / Pcij ) (11) notice that the term within the second bracket on the right-hand side of Eq. 11 is nothing but Po,ij given by Eq. 7. Thus, Eq. 11may be rewritten as follows, Tcij 2 / Pcij =(Tcij 11/12 / Pcij ) P0,ij (12) Eq. 12 can further be rearranged as follows, Tcij 2 / Pcij =(Pcij / Tcij) 7/3 P0,ij (13) By introducing Eq. 13 back into Eq. 10the following expression for Po,m is obtained, Po,m =[ΣiΣjxixj(Pcij / Tcij) 7/3 P0,ij] [ΣiΣjxixj(Tcij / Pcij)] 7/12 (14) Eq. 14 represents the general van der Waals mixing rule proposed in this work for calculation of the temperature-independent parameter, Po,m, for a multicomponent mixture. The crossparameters may be calculated using the following expression, AIChE Journal, Vol. 44, No. 10, pp. 2324-2332, 1998 7 Po,ij = (1 mij) (P0,ii.P0,jj) 1/2 (15) where mij is a binary interaction parameter. The cross-parameters for temperature and pressure in Eq. 14 may be calculated using Eq. 6b and 6d, respectively. For the case of binary mixtures, the equations to calculate the temperature-independent parameters, Po,m and P V o,m, for the liquid phase is the following, The temperature-independent parameter P v o for the vapor-phase mixture is calculated by replacing the xi’sin Eq. 16 by the compositions of the vapor phase (i.e. yi’s). The crossparameters are then given by, Po,12 = (1 m12) (P0,1.P0,2) 1/2 (16a) Tc12 = (1 k12 RM) (Tc1 Tc2) 1/2 (16b) Pc12 = 8 Tc12 / [(Tc1 / Pc1) + (Tc2 / Pc2)] (16c) Eq. 16 may be used along with Eq. 4 to predict the surface tension of binary mixtures. In the following sections the results for 55 binary mixtures are reported. Table 1 shows the physical properties of the pure fluids involved in these binary mixtures. The most widely used mixing rules for the temperature-independent parameter, Po, in Eq. 4 used in the petroleum industry were proposed by Weinaug and Katz (1943) as follows, Po,m= Σ(xi Po,i) (17) Because of the simplicity of this mixing rule, it was decided to investigate the results obtained with Eq. 4 and this simple mixing rule as a preliminary step. Thus, the surface tension of all 55 binary mixtures were predicted using the following procedure (CASE 1) (a) The temperature-independent parameter Po for the pure compounds were obtained by fitting Eq. 1 to the experimental surface tension for the pure compounds. (b) Routine VLE calculations were performed using a computer package developed in our laboratory to obtain the vaporand liquid-phase compositions as well as the saturation pressure. The results from the VLE calculations were used along with the Riazi-Mansoori AIChE Journal, Vol. 44, No. 10, pp. 2324-2332, 1998 8 equation of state to predict the vapor and liquid densities. (c) Finally, Eqs. 4 and 17 were used to predict the surface tension for the mixture. The results obtained for this CASE 1 are reported in Table 2. Note that the surface tension of all 55 sets of binary-mixture data can be predicted with an overall 2.64 AAD%. Surface Tension Prediction for Binary Mixtures As mentioned earlier, Eqs. 4, 16, 16a-16c,proposed in this work, represent the necessary expressions to predict the surface tension of binary liquid mixtures. These expressions relate the surface tension of the mixture to the bulk properties of the system such as the vaporand liquidphase compositions and vaporand liquid-phase densities. Therefore, only the temperatureindependent parameter Po for the pure components is needed. This parameter can be obtained by fitting Eq. 1 to the experimental surface tension data for the pure fluids. It may also be predicted using the corresponding-states expression given in Eq. 3. In this work, we have predicted the surface tension for 55 binary mixtures and compared the results against the experimental data for these systems. The selected systems represent a wide variety of binary mixtures (symmetric, asymmetric, and slightly asymmetric). The interaction parameter, mij, appearing in Eq. 15 has been treated as an adjustable parameter for all the mixtures investigated. The surface tension for all the binary mixtures was predicted following the procedure outlined below, (a) The temperature-independent parameter Po for the pure compounds were obtained in two ways: CASE 2 by fitting Eq. 1 to the experimental surface tension for the pure compounds; and CASE 3 by predicting them using Eq. 3. (b) VLE calculations were performed using a computer package developed according to what was discussed above to obtain the vaporand liquid-phase compositions as well as the saturation pressure. The results from the VLE calculations were used along with the RiaziMansoori equation of state to predict the equilibrium vapor and liquid densities. (c) Finally, Eqs. 4, 16, and 16a-16c were used to predict the surface tension for the mixture. The results obtained for these two cases are also reported in Table 2. It must be pointed out that the interaction parameter, mij’s in CASE 2 and CASE 3 are different. Therefore, a different value for each case is reported in Table2. Results and Discussion for Binary Mixtures Table 2 contains the results obtained for CASES1, 2, and 3. This table also reports temperature, experimental data reference, and number of data points for each binary mixture. CASE 2 The average absolute percent deviation (AAD%) obtained when the binary interaction AIChE Journal, Vol. 44, No. 10, pp. 2324-2332, 1998 9 parameter (mij) was set equal to zero is reported. Note that the surface tension for all the 55 binary mixtures can be predicted within an overall 2.06 AAD%. It may also be noticed that even when mij is set equal to zero the surface tension for most binary mixtures can be predicted with reasonable accuracy. This represents an improvement over the simple mixing rule (CASE 1)as shown in Table 2. In a different calculation for CASE 2, the binary interaction parameters, mij, were fitted to all binary mixtures. The best value for this parameter is also reported in Table 2. It may be noted that this interaction parameter is small for all the systems investigated. Note also that when this optimized interaction parameter was used the surface tension of all 55 binary mixtures can be predicted within an overall0.50 AAD%. Figures 1 through 9 show the results obtained for CASE 2 when the best value for the binary interaction parameter used in the prediction of surface tension for the binary systems. Note the excellent agreement with the experimental data. CASE 3 The AAD% obtained when the binary interaction parameter (mij) was set equal to zero and the parameters Po are predicted using Eq. 3 is reported in Table 2. Note that the surface tension for all the55 binary mixtures can be predicted within an overall 3.70 AAD%. It must be pointed out that by setting the binary interaction parameters equal to zero the surface tension is purely predicted since the only adjustable parameter is eliminated. Note that predictions thus made are in good agreement with the experimental data except in a few cases where Eq. 3 fails to yield good estimates for the parameter Po. In another set of calculation for CASE 3 the binary interaction parameter (mij) was fitted to the experimental data. The best values for this parameter are also reported in Table2. Note that that when this optimized interaction parameter was used the surface tension of all the 55 binary mixtures can be predicted within an overall 2.10 AAD%. It should be pointed out that in the absence of surface tension data for pure compounds, Eqs. 3, 4, and 14 (CASE 3) represent an excellent tool to predict the surface tension of mixtures of organic compounds. Good predictions will be obtained even when the binary interaction parameter, mij, is set equal to zero. Conclusions Based on the 55 sets of binary-mixture surface tension data analyzed in this work we may conclude that the expressions proposed in this work for surface tension prediction of mixtures are very good. These predictions are made using only the values of the temperature-independent parameter Po for the pure components in the mixture. This parameter may be obtained from experimental surface tensions of the pure components or estimated using Eq. 3. It is also worth noting that in these calculations only one adjustable parameter has been used namely the binary interaction parameter, mij, appearing in Eq. 15. There exist a number of other methods for prediction of the surface tension of liquid mixtures. None of the other methods are as accurate in predicting mixture surface tension except AIChE Journal, Vol. 44, No. 10, pp. 2324-2332, 1998 10 the ones by Carey et. al., (1980) and Winterfeld et al., (1978). The latter two methods report very accurate prediction of surface tension of mixtures provided pure fluid surface tension data are available. It may be concluded that the expressions proposed in this work are general since they have been applied to a variety of binary mixtures and the results obtained are in good agreement with the experimental data. Application of the present method to multicomponent mixtures is straight forward. However, little reliable data are available for multicomponent surface tensions for reliable comparisons. AIChE Journal, Vol. 44, No. 10, pp. 2324-2332, 1998 11 TABLE 1: Properties of the pure fluids involved in the 55 binary mixtures studied in this work _________________________________________________________________________________________ Compound Tc Pc Zc ω R* Tb ________________________________________________________________________________________ Acetone 508.1 47. 0.232 0.304 2.316 329.2 Acetonitrile 545.5 48.3 0.184 0.321 1.585 354.8 Benzene 562.2 48.9 0.271 0.212 3.748 353.2 Carbon Disulfide 552. 79. 0.276 0.109 3.079 319. Carbon Tetrachloride 556.4 49.6 0.272 0.193 3.784 349.9 Chloroform 536.4 53.7 0.293 0.218 3.071 334.3 Cyclopentane 511.7 45.1 0.275 0.196 3.310 322.4 Cyclohexane 553.5 40.7 0.273 0.212 3.966 353.8 Cis-Decalin 702.3 32. 0.245 0.286 6.281 468.9 Trans-Decalin 687.1 31.4 0.245 0.270 6.341 460.5 n-Decane 617.7 21.2 0.249 0.489 6.915 447.3 Dichloromethane 510. 63. 0.265 0.1916 2.339 313. n-Dodecane 658.2 18.2 0.240 0.575 8.269 489.5 Ethyl Acetate 523.2 38.3 0.252 0.362 3.186 350.3 Ethyl Ether 466.7 36.4 0.262 0.281 3.219 307.6 n-Hexane 507.5 30.1 0.264 0.299 4.281 341.9 Iodomethane 528. 65.9 N/A 0.1493 2.788 315.7 Iso-Octane 544. 25.7 0.266 0.303 5.619 372.4 Nitromethane 588. 63.1 0.208 0.310 1.787 374.3 Phenol 694.2 61.3 0.240 0.438 4.001 455. Tetrachloroethylene 620.2 47.6 0.250 0.254 4.341 394.4 Toluene 591.8 41. 0.263 0.263 4.450 383.8 o-Xylene 630.3 37.3 0.262 0.310 5.124 417.6 _____________________________________________________________________________________________ The critical properties and acentric factor were obtained from Reid et al., 1986. The parameter Po was obtained from experimental surface tension data for pure components; The dimensionless molar refraction R* was obtained from data for molar refraction reported by Hall (1986) or calculated using the Lorentz-Lorenz function. For Iodomethane, the acentric factor was estimated by predicting the vapor pressure at Tr = 0.7 using the Antoine equation for vapor pressure prediction. The constants for this equation were obtained from Hall (1986). For Iodomethane there was no experimental critical compressibility factor. AIChE Journal, Vol. 44, No. 10, pp. 2324-2332, 1998 12 TABLE 2 Results obtained for the 55 sets of binary mixtures investigated in this work __________ _____________________ ______________________ CASE 1 CASE 2 CASE 3 __________ _____________________ ______________________ AAD% fitted AAD% fitted Binary system AAD% mij=0 mij AAD% mij=0 mij AAD% T (K) Ref. # data AcetoneChloroform 0.42 0.48 0.0050 0.09 1.69 -0.0095 0.94 291.15 A 5 Acetone Carbon disulfide 3.06 1.42 0.0157 0.21 3.13 -0.0087 2.31 288.15 A 5 Acetone Phenol 3.64 3.64 0.0305 1.72 10.24 0.0717 3.13 293.15 A 5 Benzene Acetone 0.75 0.69 -0.0061 0.29 5.61 -0.0326 2.24 298.15 B 9 Benzene Carbon tetrachloride 0.22 0.22 0.0000 0.22 1.98 -0.0064 1.67 298.15 B 11 Benzene Carbon disulfide 2.77 0.31 0.0035 0.16 1.94 0.0000 1.94 291.15 A 11 Benzene Ethyl acetate 1.00 0.69 0.0060 0.13 6.35 -0.0346 2.79 298.15 B 11 Benzene Ethyl ether 2.46 2.36 0.0202 0.44 2.21 0.0150 1.19 291.15 A 5 Benzene n-Dodecane 1.74 6.17 -0.0705 0.70 7.06 -0.0659 0.96 298.15 C 5 Benzene n-Dodecane 2.37 4.79 -0.0705 1.72 6.16 -0.0659 1.06 303.15 C 5 Benzene n-Dodecane 1.78 4.94 -0.0705 1.30 6.01 -0.0659 1.08 308.15 C 5 Benzene n-Dodecane 2.07 6.06 -0.0705 1.05 6.88 -0.0659 0.71 313.15 C 5 Benzene n-Hexane 4.04 3.11 0.0275 0.97 2.16 0.0140 1.12 298.15 C 7 Benzene n-Hexane 3.70 2.86 0.0275 0.39 2.36 0.0140 1.04 303.15 C 5 Benzene n-Hexane 3.86 2.93 0.0275 0.69 2.38 0.0140 1.01 308.15 C 5 Benzene n-Hexane 3.18 2.39 0.0275 0.63 1.80 0.0140 1.04 313.15 C 5 Benzene Nitromethane 3.47 3.69 0.0220 1.13 2.84 0.0059 2.82 298.15 A 6 Benzene o-Xylene 0.45 0.73 -0.0060 0.05 2.40 -0.0143 0.73 298.15 D 10 Benzene Phenol 3.50 3.52 0.0338 0.28 10.67 0.0672 4.38 308.15 A 5 CCl4-Ethyl acetate 0.69 0.50 0.0020 0.45 4.87 -0.0232 2.73 298.15 B 11 Carbon disulfide – CCl4 2.55 0.76 0.0077 0.31 4.37 0.0296 1.45 293.15 E 10 Carbon disulfide CCl4 2.57 0.75 0.0077 0.25 4.39 0.0296 1.40 298.15 E 10 Carbon disulfide CCl4 2.76 0.94 0.0077 0.21 4.52 0.0296 1.43 303.15 E 10 Carbon disulfide CCl4 2.81 0.98 0.0077 0.23 4.56 0.0296 1.49 308.15 E 10 Carbon disulfide CCl4 2.75 0.92 0.0077 0.30 4.47 0.0296 1.42 313.15 E 10 Carbon disulfide CCl4 2.91 1.08 0.0077 0.41 4.59 0.0296 1.43 318.15 E 10 Carbon disulfide Benzene 3.56 1.29 0.0073 0.53 3.08 0.0065 2.95 298.15 B 6 Cyclopentane CCl4 0.29 0.31 0.0033 0.10 1.61 0.0083 0.96 298.15 D 10 Cyclopentane Benzene 2.42 2.42 0.0182 0.19 1.31 0.0070 0.77 298.15 D 11 Cyclopentane Toluene 1.88 1.27 0.0099 0.12 0.62 0.0035 0.44 298.15 D 12 Cyclopentane-Tetrachloroethylene 2.01 1.84 0.0151 0.19 2.06 0.0159 0.68 298.15 D 11 Cyclohexane Benzene 2.36 1.99 0.0169 0.14 0.82 0.0026 0.64 293.15 D 11 Cyclohexane Benzene 1.98 1.77 0.0169 0.31 1.00 0.0026 0.98 303.15 D 11 Cyclohexane cis-Decalin 0.67 1.12 -0.0102 0.12 3.70 -0.0219 1.27 298.15 D 11 Cyclohexane trans-Decalin 0.31 1.38 -0.0113 0.11 1.76 -0.0121 0.37 298.15 D 11 Cyclohexane-Tetrachloroethylene 1.95 1.94 0.0152 0.16 1.89 0.0134 0.80 298.15 D 11 Cyclohexane Toluene 1.93 1.86 0.0152 0.14 0.90 0.0055 0.61 298.15 D 13 Ethyl ether Carbon disulfide 5.27 2.49 0.0183 1.87 5.19 0.0459 4.03 291.15 A 6 Iodomethane CCl4 2.15 1.72 0.0118 0.99 2.90 0.0036 2.77 288.15 F 5 Iodomethane CCl4 2.65 2.15 0.0118 0.65 2.73 0.0036 2.60 293.15 F 8 AIChE Journal, Vol. 44, No. 10, pp. 2324-2332, 1998 13 Iodomethane CCl4 2.57 1.93 0.0118 0.51 2.34 0.0036 2.34 298.15 F 14 TABLE 2 (continued) __________ _____________________ ______________________ CASE 1 CASE 2 CASE 3 __________ _____________________ ______________________ AAD% fitted AAD% fitted Binary system AAD% mij=0 mij AAD% mij=0 mij AAD% T (K) Ref. # data Iodomethane CCl4 2.64 1.99 0.0118 0.48 2.50 0.0036 2.47 303.15 F 11 Iodomethane CCl4 2.37 1.80 0.0118 0.37 3.08 0.0036 3.21 308.15 F 9 Iso-Octane Benzene 3.55 1.18 0.0096 0.58 3.00 0.0169 1.42 303.15 G 9 Iso-Octane Cyclohexane 1.08 0.24 0.0000 0.24 2.50 0.0101 1.50 303.15 G 5 Iso-Octane n-Dodecane 1.16 1.59 -0.0147 0.29 1.53 -0.0050 1.45 303.15 G 9 Acetonitrile Carbon tetrachloride 1.56 1.97 0.0160 0.82 4.57 -0.0066 4.47 298.15 F 9 Acetonitrile Carbon tetrachloride 1.45 1.98 0.0160 0.79 4.98 -0.0066 5.02 303.15 F 8 Acetonitrile Carbon tetrachloride 1.75 2.10 0.0160 0.65 4.74 -0.0066 4.64 308.15 F 9 Acetonitrile Carbon tetrachloride 1.78 2.16 0.0160 0.54 5.59 -0.0066 5.48 313.15 F 10 Acetonitrile Carbon tetrachloride 1.86 2.22 0.0160 0.61 5.68 -0.0066 5.58 318.15 F 9 Dichloromethane-Carbon disulfide 3.21 2.97 0.0290 0.27 4.31 0.0304 3.04 293.15 H 8 Dichloromethane-Carbon disulfide 3.76 3.48 0.0290 0.20 4.06 0.0304 2.80 298.15 H 9 Dichloromethane-Carbon disulfide 3.97 3.69 0.0290 0.55 4.69 0.0304 3.15 303.15 H 10 Dichloromethane-Carbon disulfide 3.76 3.47 0.0290 0.46 4.54 0.0304 3.32 308.15 H 10 Overall AAD% 2.64 2.06 0.50 3.70 2.10 _________________________________________________________________________________________ REFERENCES: A : Meissner and Michaels (1949); B : Shipp (1970); C : Schmidt et al., (1966); D : Lam and Benson (1970); E : Luengo et al., (1988); F : Texeira et al., (1992); G : Evans and Clever, 1964; H : Aracil et al., 1989. AIChE Journal, Vol. 44, No. 10, pp. 2324-2332, 1998 14 16 17 18 19 20 21 22 23 24 25 26 27 Su rf ac e Te ns io n (d yn es /c m ) 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 x (Benzene) 16 17 18 19 20 21 22 23 24 25 26 27 Su rf ac e Te ns io n (d yn es /c m ) 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 x (Benzene) 313.15°K 308.15°K 17 18 19 20 21 22 23 24 25 26 27 28 29 Su rf ac e Te ns io n (d yn es /c m ) 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 17 18 19 20 21 22 23 24 25 26 27 28 Su rf ac e Te ns io n (d yn es /c m ) 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 x (Benzene) x (Benzene)